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ABSTRACT 


There has been much conjecture and debate about the 
possible impacts of global warming on tropical cyclone (TC) 
activity. Prior studies have presented evidence that 
global warming has already had clear impacts on TCs, that 
it has not, and that existing data sets may not be 
sufficient for determining what impacts, if any, may have 
occurred (e.g., Chan and Liu 2004; Emanuel 2005; Pielke et 
al. 2005; Webster et al. 2005; Curry et al. 2006; Klotzbach 
2006; Mann and Emanuel 2006; Wu et al. 2006; Landsea 

2007). There is even less consensus about the future 
impacts of global warming on TC activity. 

In this study, we have investigated the relationships 
between global warming and TC activity in the western North 
Pacific (WNP) region, the world's most active TC region. 
Our hypothesis was that global warming impacts on TC 
activity will occur through global warming induced changes 
in the large scale environmental factors (LSEEs) that are 
known to be important in determining the formation and 
intensity of TCs. The five LSEEs on which we focused our 
analyses were: 

Sea surface temperature (SST) exceeding 26° C 
Weak vertical shear in horizontal winds 
Large positive absolute vorticity at low levels 
Mean upward motion 
High mid-level humidity 

Our primary data sets were the JTWC TC best track 
data, NCEP reanalysis fields, and Hadley Centre SSTs. Our 
study period was 1970-2005, with a focus on late May- 
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December, the months with most TC activity in the WNP. We 
separated the data into weekly averages for 5°x5° latitude- 
longitude regions spanning the WNP. 

Using least squares linear regression, we identified 
global warming signals in both the SST and vertical wind 
shear data across the entire WNP (long term increase in 
SST, long term decrease in tropical vertical wind shear 
magnitude). These signals varied significantly on a 5°x5° 
regional scale. We also found that across the WNP, TC 
formations and accumulated cyclone energy (ACE, a measure 
of net TC intensity), underwent a long term increase during 
the study period. We then used logistic regression to 
determine to what level of significance each of the LSEFs 
affects TC formation probability. Multiple linear 
regression was performed to determine the relationship of 
each of the LSEEs to ACE. Through the two regression 
models, we determined that each of the LSEEs is significant 
for both TC formation and ACE. We then used data 
independent from the data used in the regression modeling 
to validate the models, with positive results. 

Our cumulative results support our hypothesis, and 
indicate that global warming has increased TC numbers and 
intensities in the WNP by way of the LSEFs. We will be 
conducting additional studies to determine the extent to 
which our results are sensitive to different methods of 
conducting best track analyses (cf. Wu et al. 2006). 

Additional research is also needed to better assess the 
adequacy of best track and other TC data sets for 
identifying the impacts of global warming and other long 
term climate variations on TC activity (cf. Landsea 2007) . 
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EXECUTIVE SUMMARY 


The purpose of this research is to identify the 
relationships between global warming and tropical cyclone 
(TC) activity in the Western North Pacific (WNP). The 
consensus among researchers who study TC activity is that 
TC activity (formation and strength) is a function of 
several large scale environmental factors (LSEFs) . These 
factors are: 

• Large values of low level absolute vorticity 

• Weak vertical shear of the horizontal winds 

• Sea surface temperatures (SSTs) exceeding 26 C 

• Mean upward motion and 

• High mid-level humidities 

Our study focused on three main points. We wanted to 
verify that the LSEFs given above really do influence TC 
activity. We also wished to examine whether global warming 
is exerting a detectable influence on these LSEFs, and 
subsequently, whether a global warming induced change in 
these LSEFs is exerting a detectable influence on TC 
activity. 

Data for these factors was gathered for the years 1970 
(the beginning of the satellite era, ensuring our TC count 
is reasonably correct and consistent) through 2005 and 
averaged into a weekly grid of five by five degree blocks 
covering the WNP. Regression was performed to determine 
how these factors are changing with time. We determined 
that SSTs are increasing, as is shear, though in the 
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tropics, where TC activity is greatest, the magnitude of 
shear is actually decreasing. 

Logistic regression was used to examine the 
relationship between the LSEFs and the probability of TC 
formation. It was determined that all the previously 

mentioned LSEFs are important in influencing TC formation 
probability, with SST and vertical velocity being the 

factors to which TC formation is most sensitive. Using the 
logistic regression fit based on our 1980 - 2005 data, we 
estimated the individual TC formation probabilities by week 
for each of the blocks of our grid. We then displayed 

weekly estimated probability contour plots showing the 
probability of formation with actual TC formations 

superimposed. The performance of these probability contour 
plots against actual TC formations validates the model. 
The logistic regression model for the probability of TC 
formation is: 

log(p/(l-p)) =-22.3677 + 0.0066SST + 1.35Coriolis 
+ 0.0407AbsVort + 0.1436VertVel 
-0.027GRelHum - 0.0264Shear - 0.0031Shear^ 

where p is the estimated probability of TC formation. 

Likewise, linear regression and analysis of variance 
were used to determine what if any relationship exists 
between the LSEFs and accumulated cyclonic energy (ACE, a 
measure of the strength of a storm). Again, it was 
determined that all the previously mentioned LSEFs appear 
to influence ACE, with SST and absolute vorticity being the 
factors to which ACE is most sensitive. Our model was 
validated by comparing our predicted values of ACE to 
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actual values; there was a surprisingly good agreement. 
Our validated model is: 

ACE = -92365 + 26.52SST + 3005.25AbsVort 
+ 58.83VertVel + 211.99RelHum 
+ 55.32Shear - 5.09Shear^ 

with a standard error of 20380 kts^ on 4448 degrees of 
freedom. 

Therefore it is our contention that through the LSEFs, 
global warming has and will continue to increase the 
average number and intensity of TCs in the WNP. 
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I. INTRODUCTION 


A. BACKGROUND 

Global warming is a controversial topic. Is it for 
real? Is much being made of nothing? Is combating it 

worth the required significant societal and economic 
changes ? 

The possible manifestations of global warming are 
especially controversial. Particularly after hurricane 
Katrina hit the U.S. in August 2005, there has been much 
conjecture that recent changes in hurricane activity in the 
Atlantic^ have been due to global warming^. The purpose of 
this study is to identify what (if any) environmental 
changes have been occurring that may have impacted tropical 
cyclone (TC) activity. 

Changes in TC activity (we define TC activity as TC 
formation and TC intensity) have significant ramifications 

for U.S. military basing/staging; the health of the 
economy, and the U.S. population as a whole; and there are 
certainly ramifications for the entire "global village". 

Prior studies show conflicting results. Klotzbach 
maintained there is no global warming - TC activity link^. 
Emanuel asserted that there is no discernable trend in TC 

formations per year, but that intensity is increasing^. Chan 
and Liu showed a relationship between TC intensity in the 
Western North Pacific and El Nino, but no relationship 

^ P. J. Webster, G. J. Holland, J. C. Curry, H.-R. Chang, Science 
309 , 1844 (2005). 

2 R. A. Anthes, R. W. Corell, G. Holland, J. W. Hurrell, M. C. 
MacCracken, K. E. Trenberth, Bull.Amer.Meteor.Soc 87 , 623 (2006). 

^ P. J. Klotzbach, Geophysical Research Letters, 33 , L10805 (2006) . 

^ K. Emanuel, Nature, 436 , 686 (2005) . 


1 



between intensity and global warming^. Webster et al. 
concluded that TC formations and intensity are both 
increasing®. Pielke et al. make the point that the impact 
of TC activity on mankind is dominated by societal 
vulnerability, and not by increases in TC activity'^. 
Certainly the focus of these studies has been on TC 
activity in the Atlantic. Wu et al. identified notable 
differences between several TC data sets for the WNP in 
terms of their depiction of changes in TC activity over the 
last several decades®. Landsea concluded that the available 
TC data is insufficient to allow any impacts of global 
warming on TC activity to be distinguished from variations 
in TC activity due to other causes® (e.g., natural multi- 
decadal climate variations). Trenberth and Shea asserted 
that global warming has resulted in increasing global 
average sea surface temperatures^®. 

A common characteristic of these prior studies is that 
they focused on very large spatial and temporal scales. 
Emanuel (2005), for example, calculated a basin wide 
measure of TC intensity and compared it to tropical SST 
averaged over a very large portion of the same basin. Both 
quantities were smoothed twice before correlating them. 
Our concern is that this large scale averaging and 


® J. C. L. Chan, K. S. Liu, Journal of Climate, 17 , 4590 (2004) . 

® P. J. Webster, G. J. Holland, J. A. Curry, H.-R. Chang, Science, 
309 , 1844 (2005). 

® R. A. Pielke Jr, C. Landsea, M. Mayfield, J. Laver, R. Pasch, 
Bull.Amer.Meteor.Soc 86 , 1571-1575 (2005) 

® M.C. Wu, K.H. Yeung, W.L. Chang, Eos, 87 , 537-538 (2006) . 

® C.W. Landsea, Eos, 88 , 197-208 (2007) . 

K. E. Trenberth, Dennis J. Shea, Geophysical Research Letters, 
33 , (2006). 
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smoothing significantly dilutes the information necessary 
to understand changes in TC activity related to global 
warming. 

1. Influences on Tropical Cyclone Activity 

Over the past 50 years, significant strides have been 
made in the understanding of tropical cyclone activity. The 
consensus is that TC formations, intensities and tracks, 
are strongly influenced by several primary large scale 
environmental factors (LSEF) , includingi^: 

• Sea surface temperatures (SSTs) exceeding 26 C 

• Weak vertical shear of the horizontal winds 

• Large positive absolute vorticity at low levels 

• Mean upward motion 

• High mid-level humidity 

The LSEFs are used to analyze and forecast TC 
activity, especially TC formations. Some seasonal 

predictions of TC activity in the North Atlantic are based 
on derived empirical relationships, including relationships 
between the LSEFs and TC formations^^. 

In the context of trying to show a relationship 
between TC activity and global warming, the state of the 
art is to use these factors on a very grand scale. In some 
cases, SSTs, by far the most tangible of the LSEFs, are 
used almost exclusively for analyzing and predicting TC 
activity, with the remaining LSEFs being virtually ignored. 


W. M. Frank, 1987, Tropical Cyclone Formation. Chap. 3, A Global 
View of Tropical Cyclones. Office of Naval Research, Arlington, VA 
22217, 53-90. 

W. M. Gray, 1975: Tropical Cyclone Genesis in the Western North 
Pacific, Naval Air Systems Command, Washington DC, 20361. 
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2. Using the Large Scale Environmental Factors on an 
Appropriate Scale 

We think that much is lost in the understanding of TC 
activity if, as in prior studies, only one of the LSEFs 
(e.g., SST) is used to assess the impacts of global warming 
on TC activity, and if an LSEF is analyzed using data 
averaged over millions of square kilometers and several 
months, seasons, or years. The prior studies based on 
these sorts of analyses have had mixed results^^'. 
Contrary to the approach of these prior studies, we contend 
that all of the LSEFs are important in determining the 
impacts of global warming on TC activity, and that they do 
so on relatively local spatial and temporal scales (with 
respect to TC scales). In other words we think a tropical 
cyclone forms because the LSEFs on scales of several 
hundred kilometers, and over a few days or up to a couple 
of weeks are favorable for TC formation, not because the 
average temperature for the tropical North Atlantic for the 
year 2002, smoothed out over five years exceeded 26° C. 

Therefore, we decided to study the relationship 
between global warming and each of the LSEFs to determine 
if any predictions regarding TC activity can be made given 
the effect of global warming on each LSEF, if indeed there 
appears that there is such an effect. We also decided to 
analyze the LSEF and TC data at space and time scales 
consistent with the scales at which the LSEFs are thought 
to influence TC activity (i.e., scales of about a week and 
hundreds of kilometers, rather than scales of a whole TC 
season (about six months) and a whole tropical basin 
(thousands of kilometers) ) . Thus in this study, we test 

K. Emmanuel, Nature, 436 , 686-688 2005). 

Johnny C. L. Chan, Science, 311 2006). 
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the hypothesis that if global warming is influencing TC 
activity, it is doing so through all the LSEFs, and it is 
doing so at scales of hundreds of kilometers and about a 
week. 

Most of the data on the LSEFs is available in six hour 
time intervals and in two and a half by two and a half 
degree latitude by longitude geographical increments (about 
250 by 250 km squares) . This allowed us to analyze the 
LSEFs on scales that fit their known impacts to TC 
activity. This should allow us to confirm that the LSEFs 
have an impact on TC activity, if they do. 

The main purpose of the regression models is to 
facilitate the identification of relationships between 
global warming and TC activity, through the LSEFs commonly 
understood to play a role in TC activity (Figure 1). 
However, it is possible that these models might also be 
useful in analyzing and forecasting TC activity in general. 
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Figure 1. Proposed regression modeling. 

B. SCOPE, LIMITATIONS, AND ASSUMPTIONS 

1. Scope of the Work 

This thesis will build on the prior studies of the 
climatology of TC activity, and the relatively limited work 
done investigating TC-global warming relationships. 

The research conducted will focus on the Western North 
Pacific (WNP) for many reasons. The WNP is: 

1. the region with the largest annual number of TCs 
(approximately 35 percent of the global total) 

2. a region where TC activity has demonstrated 
effects (via teleconnections) on the U.S., and 
even Europe 

3. bordered by nations which comprise a large 
portion of the Earth's population 

4. a region of interest to the U.S. military 


6 










Data used will include all available data during 1970- 
2006, from the onset of extensive satellite coverage to the 
present. Data prior to 1970 will not be used, since prior 
to full satellite coverage, there were many likely TCs that 


never made 

landfall 

and 

for which 

little or 

no 

data 

collected. 

Before 

the 

satellite 

era, TCs 

at 

sea 


documented as they were observed by ship or aircraft, or if 
they made landfall. Between 1945 and 1970 an increase in 
reported TC formations is quite possibly a result of an 
increasing number of ships and aircraft at sea. As there 
have been many contentious arguments over the TC data base, 
we deliberately chose to use data from the satellite era 
only to avoid this issue as best as we could^^. 

2. Limitations 

There are two choices for measures of global warming 
to be used as independent variables, atmospheric CO 2 
concentration in parts per million (ppm) (see Figure 2), or 
globally averaged surface temperature anomaly (GASTA). CO 2 
concentration is a rather direct cause of global warming, 
but it does not capture the whole story of greenhouse 
emissions. GASTA of course is not a cause of global 
warming, but rather an attempt to quantify it. The 
limitation of GASTA is that while there is little debate as 
to whether Earth is experiencing a warming period, there is 
much more debate as to whether that warming is due to 
anthropogenic (man-made) causes. Figure 3 shows the 
existing relationship between CO 2 concentration and GASTA 
from 1970 through 2006. 


C.W. Landsea, Eos, 88 , 197-208 (2007) . 
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Figure 2. Time series plot of CO 2 concentration's (ppm) for 
the time frame 1970-2006, with twelve monthly data points 
per year. 


The basic science behind greenhouse gasses and global 
warming is well established and understood. By altering 
the make up of the Earth's atmosphere by adding combustion 
gasses (primarily CO 2 ) from industry, power plants, and 
vehicles, the radiation of energy from Earth to space is 
changed in a way that causes Earth to get warmer. What is 
far less understood is what happens to Earth as it gets 
warmer. Earth is an enormous mechanical system, with a 
myriad of feedback cycles that operate on scales of minutes 
to millennia, making it very difficult to predict the 
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results of adding greenhouse gasses to the atmosphere (and 
also making simplistic predictions of the results of global 
warming laughably wrong). 



Figure 3. Global Average Surface Temperature Anomaly^'^ 
(deg C multiplied by 100) vs. CO 2 Concentration (PPM). 


Figures 2 and 3 show that CO 2 concentration is 
increasing approximately linearly with time, as is GASTA, 
and that there is a linear relationship between CO 2 
concentration and GASTA. In our study, we have not 
attempted to determine the causes of global warming. Our 
work is focused on the relationships between observed 
global warming and observed TC activity in the WNP, not on 
the causes of global warming, natural or otherwise. One 
other clarification: as we have shown a linear relationship 
between CO 2 concentration and year, and a linear 

Global Average Surface Temperature Anomalies, 
http://Iwf.ncdc.noaa.gov/oa/climate/research/anomalies/anomalies.html. 
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relationship between CO 2 concentration and GASTA, we will 
use years as a proxy measure of global warming. We will 
use years explicitly as an independent variable to show its 
relationship with each of the LSEFs and with WNP region 
wide TC activity. The purpose of using years as a proxy is 
to aid the reader, since GASTA is a less tangible 
independent variable than years. 

Our study will not establish causal relationships. 
However, we do expect the evidence that we develop in this 
study will contribute to a determination of the causes of 
any long term changes in the LSEFs and in TC activity in 
the WNP. The determination of these causes will require 
the application of many different methods, including 
statistical analyses of observational data (such as ours), 
theoretical studies, and modeling studies. 

Our study will focus on identifying the relationships 
between the LSEFs and (a) our proxy measure of global 
warming and (b) TC activity. For our analysis of the 
effects of the LSEFs on TC activity, the LSEFs will be our 
direct independent variables and TC formation and 
accumulated cyclone energy (ACE) will be our response 
variables. ACE quantifies TC intensity by squaring TC wind 
velocity and summing at six hourly intervals. For our 
analysis of the impacts of global warming on TC activity, 
the LSEFs will serve as our intermediate response 
variables, with our final response variables being TC 
frequency and ACE. It is important to keep in mind that 
the statistical relationships we will show do not by 
themselves allow us to determine causality. However, the 
relationships may contribute to the determination of 
causality, as discussed in the preceding paragraph. 


10 



An example is in order. A regression analysis between 
atmospheric CO 2 concentration (which has been steadily 
increasing during 1970-2006) and GASTA over the time frame 
(using the statistical software package S plus) gives: 

Temp(°C*100) = -358.25 + 1.091C02Concentration (1) 
with a standard error of 12.28 on 1922 degrees of freedom. 

This result estimates that for every 1 ppm increase in 
atmospheric CO 2 , mean GASTA increases by 0.011 degrees C, 
with a standard error of 0.017. The adjusted is 0.69 
which measures the strength between GASTA and CO 2 
concentration. 

The result is a surprisingly strong relationship that 
might suggest irrefutable proof that greenhouse gasses are 
causing an undeniable warming of the planet; certainly the 
strength of the relationship ought to be part of the global 
warming discussion. Figure 4 however, shows the 
relationship between the author's age (which unfortunately 
also has been steadily increasing with time) and GASTA. 
Regression analysis exploring the relationship between the 
author's age and GASTA yields the following equation: 

Temp(°C*100) = -15.5 + 1.68Age (2) 
with a standard error of 12.6 on 1922 degrees of freedom. 

Now the results estimate that for every year increase 
in my age, mean GASTA increases by 0.017 degrees C, with a 
standard error of 0.026 and an adjusted R^ of 0.67. While 
the relationship is strong, it is unlikely that my age is 
causing global warming. 
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Figure 4. GASTA vs. the author's age (years). 

In other words, regression can show the strength of a 
relationship between predictor and response variable, but 
additional information is needed to establish whether the 
relationship is physically plausible and therefore 
indicative of causality. In this case, reasoning based on 
additional information from atmospheric science theory and 
modeling indicates that global warming over several decades 
may influence TC activity. Therefore, a statistical 
relationship between global warming and TC activity may be 
physically plausible and indicative of causality; thus the 
relationship is relevant to the determination of causality. 
However, there are no theoretical or modeling indications 
that the author's age is related to changes in TC activity. 
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We will spare the reader the regression we did showing an 
inverse relationship between hair follicle density on his 
head and global warming. 

3. Assumptions 

Initial assumptions and hypotheses made in support of 
this work include the following: 

a. Assumptions 

1. That the spatial and temporal resolution of the 
data is sufficient to answer our questions. 

2. Our global warming index, years, sufficiently 
captures the major signals of global warming, 
and impacts (or does not impact) the LSEFs, 
which in turn influence TC activity, in such a 
way as to be useful when using regression 
tools. 

Our work will extend the prior research with the goal 
of answering the following questions with analysis based on 
local values of the LSEFs and TC activity rather than 
values averaged over large regions of both time and space 
as has been done in previous work. 

1. In the WNP, does there appear to be a 

relationship between global warming and any of 
the LSEFs which influence TC activity? 

2. In the WNP, does there appear to be a 

relationship between global warming and TC 
activity? If so, for what aspects does global 
warming appear to have impacts? 

3. Can a relationship between the LSEFs be modeled 
and used to establish relationships between 
global warming and TC activity? 

4. If relationships between global warming and the 

LSEFs can be identified, and if the 

relationships between the LSEFs and TC activity 
can be adequately modeled, then can the impacts 
of global warming on TC activity be modeled and 
potentially predicted? 

b . Hypotheses 

The above questions will be addressed in the 
process of exploring the following hypotheses: 
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1. Any relationship between global warming and TC 
activity is because global warming is 
influencing the LSEFs, and because the LSEFs 
influence TC activity. 

2. TC activity is a function of local values of 
the LSEFs not values averaged over large 
regions of time or space. 

3. TC formations are a function of each of the 
LSEFs. 

4. ACE is a function of the each of LSEFs. 

Hypothesis testing will be performed as follows: 

1. Using least squares we will show the nature of 
the relationship between each of the LSEFs 
(each being a dependent variable) and global 
warming (where year is our proxy and the 
independent variable). This will support the 
hypothesis that global warming influences the 
LSEFs. 

2. Using the tools of regression, a TC formation 
probability model will be developed using local 
values of each of the LSEFs. The model will be 
validated by comparing performance of the model 
to actual TC formation locations. 

3. Using the tools of regression, an ACE model 

which estimates ACE given a TC has formed will 
be developed using local values of each of the 
LSEFs. The model will be validated by 

comparing ACE totals estimated by the model to 
actual ACE totals for each year. 

C. OUTLINE 

Chapter II briefly discusses the data sources 
available and the data sets used to complete this work, as 
well as the rationale as to why particular data sources 
were chosen. The methods for error checking the raw data 
and converting the data into a data set designed for ease 
of using regression techniques are also discussed. Sample 
plots to demonstrate error checking are included. Finally, 
we define the WNP region for the purposes of our work, and 
explain the rationale for choosing this area. 
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Chapter III begins with addressing what large scale 
region wide changes can be seen in the LSEFs. These 
changes are then revisited on a smaller geographical scale 
to determine if the observed changes are uniform or have 
some sort of regional dependence. Large scale region wide 
trends in TC activity, both frequency of formation and ACE, 
also are addressed. Eollowing the large scale region wide 
analyses, the results from regression models for TC 
formation probability and ACE are discussed, along with 
model validation results. Simple sensitivity analysis 
results for TC formations and ACE are shown as well. 

In Chapter IV we outline the many conclusions and 
recommendations we developed from the model results 
presented in Chapter III. Areas for further research are 
discussed as well. 
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II. DATA AND METHODS 


A. STUDY PERIOD AND REGION 

The time period we selected for our analysis is 1970 
through 2006. This time period was chosen to maximize the 
size of the data set while utilizing data only from the WNP 
satellite era. While TCs form year round in the WNP, 
formations peak in the months of July through October and 
are at their minimum in Februaryi^ so to control the size of 
our data set we only used data from weeks 20 through 52 
(mid-May through December). 

The region selected for our study was 0-40 North and 
105-185 East. This region was chosen because it is 
consistent with areas defined as WNP by other researchers, 
it minimizes the area that could include Eastern North 
Pacific TCs, and it allows for the common northward 

recurvature of many WNP TCs. 

B. DATA SOURCES AND SELECTION 

There is a large amount of data available to quantify 
both the LSEEs, and TC activity. Table 1 shows a listing 
of data sources, with the data type, periodicity, temporal 
coverage, and spatial coverage. Table 2 shows data used in 
a selection of prior studies of the links between global 
warming and TC activity. 


J. Camargo, A. W. Robertson, S. J. Gaffney, P. Smyth, 

WWW.datalab.uci.edu/papers/camargo_etal.pdf 

J. C. L. Chan, K. S. Liu, Journal of Climate, 17 , 4590 (2004). 
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Data set 

Coverage 

(years) 

Grid 

Resolution 

SST 

Velocity 

Vorticity 

Humidity 

Upward 

Velocity 

NCEP I 

1948-2006 

Lat/lon 

2.5X2.5 deg, 

6 hourly 


X 

X 

X 

X 

NCEP II 

1979-2006 

Lat/lon 

2.5X2.5 deg, 

6 hourly 


X 

X 

X 

X 

(NCEP) 

OISST 

1982-2006 

Lat/lon 

1X1 deg, 6 
hourly 

X 





(NCEP) 

ERSST 

1854-2006 

Lat/lon 

2x2 deg, 
monthly 

X 





Hadley 

Centre 

1854-2006 

Lat/lon 

1x1 deg, 
monthly 

X 






Table 1. A selection of potential LSEF data sources. 


Prior studies 

LSEF data 

TC activity data 

Chan and Liu (2004) 

NCEP 

JTWC 

Webster et al. 

(2005) 

Not cited 

Not cited 

Emanuel (2005) 

Hadley Centre 

JTWC 

Trenberth and Shea 
(2006) 

Hadley Centre 

N/A (Atlantic 

focus) 

Klotzbach (2006) 

NCEP 

JTWC 

Mann and Emanuel 

(2006) 

Hadley Centre 

N/A (Atlantic 

focus) 


Table 2. A selection of prior studies and their cited data 

sources. 


Part of our approach in selecting the data for our 
study was to use, as much as possible, data that had been 
used in the key prior studies of global warming links to TC 
activity, in order to allow relatively direct comparisons 
of our methods and results to these prior studies. 

We chose two measures for global warming, CO 2 
concentration (see Figure 2) obtained from the National 
Oceanic and Atmospheric Administration^o, and global average 


Trends in Atmospheric Carbon Dioxide - Mauna Loa, 
http://WWW.cmdl.noaa.gov/ccgg/trends/. 
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surface temperature anomaly (GASTA), obtained from the 
National Climatic Data Center^i. Figure 5 shows a time 
history of GASTA. 


(ilobai Mean I emperature over Land & ucean 



with respect to 1901-2000 mean (from NOAA). 

CO 2 concentration as measured at Mauna Loa, Hawaii, is 
available at monthly intervals and certainly is a 
potential, although simplistic, predictor for global 
warming. GASTA is also available in monthly increments, 
and while it is a direct measure of the warming of the 
Earth, the extent to which it represents anthropogenic 
forcing is unclear22. GASTA may also cause complications in 
that when used as an independent variable, it may not be 
sufficiently independent of one of our response variables, 
SST. However, based on the relationships shown in Chapter 
I, years will be used as a proxy measure of global as we 
are more interested in changes in TC activity per year than 
changes in TC activity per ppm CO 2 . 

Global Surface Temperature Anomalies, 
http://WWW.ncdc.noaa.gov/oa/climate/research/anomalies/anomalies.html. 

'^'^Intergovernmental Panel on Climate Change, Working Group I, 
http://ipcc-wgl.near.edu/wgl/wgl-report.html. 
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SST data are from the Hadley Center "ISSTl" data set^^ 


from the UK Met Office website at 
http://hadobs.metoffice.com/hadisst/. This data is at a 
spatial resolution of one degree by one degree, and a 
monthly temporal resolution. The National Centers for 
Atmospheric Prediction (NCEP)/National Center for 
Atmospheric Research (NCAR) reanalysis SST data set was not 
selected because we wanted to be consistent with prior 
studies, the Hadley Centre data was at a finer spatial 
resolution, and because the NCEP OISST data set, while at a 
one degree and six hourly resolution, was not available for 
the full study period. The Hadley Centre ISSTl data set 
has been used by other researchers (Table 2), and is at the 
same temporal and better spatial resolution as the NCEP 
ERSST data so we chose to use the ISSTl set. Eigure 6 is a 
sample contour plot of data. 


N. A. Rayner, D. E. Parker, E. B. Horton, C. K. Folland, L. V. 
Alexander, D. P. Rowell, E. C. Kent, A. Kaplan, Global analyses of 
sea surface temperature, sea ice, and night marine air temperature 
since the late nineteenth century J. Geophys. Res.Vol. 108, No. D14, 
4407 10.1029/2002JD002670. 
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Figure 6. August 1997 SST contour plot (degrees C*100). 

For the atmospheric LSEFs in the WNP, we chose the 
NCEP/NCAR reanalysis fields. Two sets of reanalysis 

fields were used: Reanalysis I, which is populated by 
variables taken at 2.5 by 2.5 degree increments on six hour 
intervals starting from 1948 and continuing to present,^4 
and the Reanalysis II fields which start in 1979 and 
continue to present using the same temporal and spatial 
coverage as the Reanalysis I fields. Because the 

Reanalysis II fields are an improvement over Reanalysis 1^5, 
the Reanalysis II fields were used over the entire 
timeframe for which they are available. The Reanalysis I 
field data were only used for 1979 and earlier. Since 
there is overlap between the two data sets, we compared 

24 Kalnay et al.,The NCEP/NCAR 40-year reanalysis project. Bull. 

Amer. Meteor. Soc., 77, 437-470, 1996. 

25 NCEP-DOE Reanalysis II: Summary, 

http://WWW.cdc.noaa.gov/cdc/data.ncep.reanalysis2.html#detail. 
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data from both fields (e.g.. Figure 7) and decided that 
little if any significant error would be introduced by 
using the two different data sources for most aspects of 
our study. The NCEP Reanalysis data was provided by the 
NOAA/OAR/ESRL PSD, Boulder, Colorado, USA, from its web 
site at http://www.cdc.noaa.gov/ . 
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Eigure 7. Contour plots of zonal wind speed at 200mb level 
based on data from the Reanalysis-II (top) and Reanalysis-I 
(bottom) data sets. 
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For our measures of TC activity, we chose to use 
tropical cyclone best track data (position, date and time, 
and wind velocity)^6 from the Joint Typhoon Warning Center 
(JTWC) . The TCs included are all "wp" labeled best tracks 
from the JTWC, not just those that originate inside of our 
WNP region. There are several sources of TC data for the 
WNP, and some potentially important differences between 
these sources^^. This has led to disagreement about what 
has actually happened with TC activity over the last 
several decades^s. Our selection of the JTWC best track 
data set was based on it having been used in the 
preponderance of the literature reviewed. Figure 8 shows 
the JTWC analysis of the WNP TC tracks for 1997. 


Joint Typhoon Warning Center Products 
https://metocph.nmci.navy.mil/jtwc.php. 

M.C. Wu, K.H. Yeung, W.L. Chang, Eos, 87 , 537-538 (2006) . 

C. W. Landsea, Eos, 88 , 197-208 (2007) . 
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TC Best Tracks for year 1997 



Figure 8. 1997 TC tracks in the WNP. Each dot marks the 
chronologically ordered position of a TC in six hourly 
intervals according to the best track data set. Note the 
number of recurving TCs (those that have a strong northward 
component). Because of recurving TCs, our WNP region 

extends to 40N. 

C. VARIABLES 

1. First Stage Independent Variables for Global 
Warming/Environmental Factors 

• CO 2 concentration (ppm) 

• Global averaged surface temperature (°C) 

• Year (used as a proxy for the variables 
above based on strong linear relationships) 

2. Response Variables for First Stage Independent 
variables (year) and Intermediate Stage 
Independent Variables 

• Vorticity (s~^*10^) measured at 850 mb 

• Vertical wind shear (m/s), the 200 mb wind 
velocity minus the 850 mb wind velocity 

• Sea Surface Temperature (°C*100) 
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• Relative Humidity (%) at 500 mb 

• Vertical Velocity (-l*Pa/s) at 500mb Note 

that as air ascends, it experiences lower 
and lower atmospheric pressure, thus 
pressure (in Pascals) decreases as ascent 
occurs. We multiplied the NCEP reanalysis 

vertical velocities by -1. Thus, in our 

results, positive vertical velocities 
correspond to upward air motions. 

3. Response Variables for Second Stage Independent 

Variables 

• TC formation (number) 

• Accumulated cyclone energy (ACE, kts^) a 

measure of the maximum sustained TC winds^^. 
Only data for TCs that reached tropical 
storm intensity (maximum sustained winds of 
18 m/s (35 kts) or greater) were used in our 
calculations of ACE. 

D. TEMPORAL-SPATIAL DATA BLOCKS 

Data was obtained for the period 1970-2006 (however, 
as of the end of this study, the 2006 JTWC best track data 
was not available) . We averaged the six hourly NCEP 
reanalysis data into one week intervals. The purpose of 
this averaging was to bin the data into periods 
approximately consistent with those associated with the 
effects of the LSEEs on TC activity^^, keep the data set to 
a manageable size, and keep the program run time (for 
manipulating the data set) to a practical amount. We felt 
that operating on a weekly time scale would not unduly 
dampen the LSEE effects critical for TC activity and would 
be a vast improvement in temporal resolution over that in 
prior studies^i'32 ^ 

29philip J. Klotzbach, Geophysical Research Letters, 33 , (2006). 

3*3 w. M. Gray, 1975: Tropical Cyclone Genesis in the Western North 
Pacific, Naval Air Systems Command, Washington DC, 20361. 

33 M. A. Mann, K. A. Emanuel, Eos, 87 , 233 (2006) . 

32 j. c. L. Chan, K. S. Liu, Journal of Climate, 17 , 4590 (2004). 
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We also averaged the data into five degree latitude by 
five degree longitude blocks (about 500km by 500km) to 
represent the spatial scale at which the LSEFs are thought 
to influence TC activity^^. Thus, the data was analyzed in 
temporally and spatially averaged blocks covering one week 
and five by five degrees per block. The end result is that 
over our defined area of the WNP (see Figure 9) we have 
145,000 data blocks that include over 1000 TC formations 
and 4500 blocks where there was TC activity i.e., either a 
TC formation and/or positive ACE. 



Figure 9. TC formation sites for 1970 - 2005, with our WNP 
region shown by the red box. 


W. M. Gray, 1975: Tropical Cyclone Genesis in the Western North 
Pacific, Naval Air Systems Command, Washington DC, 20361. 
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E. STATISTICAL ANALYSES METHODS 

Several statistical methods were used in analyzing the 
data identified above. Least squares fitted lines^^ were 
used to identify any trends on a broad scale across the 
WNP. This was our tool chosen to evaluate the relationship 
between global warming and the LSEFs, and to estimate LSEF 
changes using the independent variable, years. Likewise, 
least squares were used on a broad scale to estimate the 

change in TC activity over our time period, again using 
years as our independent variable. 

In addition, trellis plots were used to re-examine the 
behavior of the LSEFs. The trellis plots are used to 
explore the same data used for the least squares fits, 
except that the data was partitioned into smaller 

geographical sub regions as opposed to being analyzed for 
the overall WNP region. By use of trellis plots, we were 
able to examine whether any changes to the LSEFs were 
occurring uniformly across the WNP or whether there was, as 
we suspected, some sort of regional dependence on how the 

LSEFs changed. Seeing a regional dependence would provide 
some indication that analyzing our data on the small, 
weekly five by five degree scale might be more sensitive 
than using larger areas and periods, as done in prior 

studies of the effects of global warming on TC activity 
(see Chapter I) . 

Logistic regression's was used to estimate the 
relationship between the log-odds of TC formation, and the 

J. L. Devore, 2003: Probability and Statistics for Engineering and 
the Sciences 6^^ ed. , Brooks/Cole-Thomson Learning, 10 David Dr., 

Belmont, CA 94002. 

D. Collett, 1991, Modelling Binary Data. Chap. 5, Overdispersion. 
Chapman and Hall, 2-6 Boundary Row, London SEl 8HN, UK, 188. 
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LSEFs. Our blocks of LSEF data were used to estimate an 
equation whereby the log-odds of TC formation probability 
is a linear function of the LSEFs. Using independent data, 
in weekly five by five degree blocks, we could then check 
and validate the performance of this model against expected 
behavior (e.g., the monsoon trough, a region in the 
tropical WNP where the LSEFs are often favorable for TC 
formation's) and against actual TC formations. 

Multiple linear regressions^ was used to estimate the 
relationship between expected ACE and each of the LSEFs. 
Our blocks of LSEF data were used to estimate an equation 
whereby the expected amount of ACE generated in a weekly 
five by five degree block given a TC has formed is a linear 
function of the LSEFs. Using independent data in weekly 
five by five degree blocks, we could then check and 
validate the performance of the model against actual ACE 
generated in a given year. These analyses result in the 
following chain: 

1. We used least squares to identify any broad scale 
relationships between global warming and the LSEFs. 
This analysis helped us to test the idea that global 
warming, if it has impacted TC activity, has done so 
via the LSEFs. 

2. We examined whether broad scale changes in the LSEFs 
had any regional dependence. This helped us test 
the idea that the impacts of global warming need to 
be analyzed at scales comparable to those at which 
the LSEFs affect TC activity. 

3. We used logistic regression applied to our small 
scale weekly five degree by five degree blocks to 


L. Elsberry, W. M. Frank, G. J. Holland, J. D. Jarrell, R. L. 
Southern, 1971: "A Global View of Tropical Cyclones", Office of Naval 
Research, Arlington, VA 22217. 

J. L. Devore, 2003: Probability and Statistics for Engineering and 
the Sciences 6^^ ed. , Brooks/Cole-Thomson Learning, 10 David Dr., 

Belmont, CA 94002. 
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estimate the relationship between the LSEFs and TC 
formation. This analysis helped us to test the idea 
that global warming, if it has impacted TC 
formations, has done so via the LSEFs. 

4. We used multiple linear regression applied to our 
small scale weekly five degree by five degree blocks 
to estimate the relationship between the LSEFs and 
ACE. This analysis helped us to test the idea that 
global warming, if it has impacted TC intensity, has 
done so via the LSEFs. 
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III. ANALYSIS AND RESULTS 


A. LONG TERM TRENDS IN THE LARGE SCALE ENVIRONMENTAL 

FACTORS 

1. Sea Surface Temperature 

In concert with the trend observed in Figure 5 and as 
discussed by Trenberth and Shears, SSTs in the WNP have 
steadily increased with time during 1970-2006(see Figure 
10) . The trend line in the graph is a least squares fit, 
with the equation for SSTs being: 

SST (°C*100) = -338.5859 + 1.5138*year (3) 

yielding an estimated 0.529 degree centigrade increase in 
mean SST during 1970-2006. The standard error for the 
intercept and the year coefficients are 184.80 and 0.093 
respectively. The coefficient for our year term, +1.5138 
is an estimate of the relationship between year and SST, 
and has units of °C*100/year. As year is our proxy for 
global warming, equation (3) indicates that there is a 
positive relationship between global warming and WNP SSTs 
(but, as discussed in Chapter I, this relationship does not 
indicate causality, although it may contribute to a 
determination of causality) . Note too, that the 

coefficient scales the equation (our SSTs are in units of 
°C*100) and keeps the units consistent, so the size of the 
coefficient is not indicative of the importance of the term 
that it scales. 


K. E. Trenberth, Dennis J. Shea, Geophysical Research Letters, 
33 , (2006). 
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Figure 10. SSTs (°C*100) averaged weekly and over all 

five degree blocks inside the WNP vs. year. Note the 
increasing trend. 


Realize that Equation 3 quantifies SST over the entire 
WNP. While Figure 10 aids in an overall understanding of 
what is happening with the SST LSEF, we remind the reader 
that we are more interested in what is happening in the 
smaller scale of weekly time frames and five degree blocks, 
since we expect that the LSEFs on these smaller scales will 
tell us more about TC activity. 

Figure 11 below shows box plots of the same overall 
WNP SST data used above. However, Figure 12 shows that the 
behavior of SST clearly is dependant on both latitude and 
longitude (week is still being ignored here for the sake of 
simplicity). Positive trends are still apparent in all 


nine regions 

of 

the 

WNP shown in Figure 

12. But the 

behavior of 

the 

SSTs 

even in these large 

regions varies 
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considerably from region to region. This indicates that 
characterizing the behavior of SST and other LSEFs using 
WNP wide averages (averages over millions of square km) can 
mischaracterize what is happening in the vicinity of TC 
formation sites or TC tracks. 



Figure 11. Box plots of SSTs over entire WNP using data 

for weeks 20-52. 
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Figure 12. Conditional plot showing SSTs conditioned on 

latitude and longitude. Increasing trend in SST over time 
is evident, though different from region to region. The 
temperatures marked with green backgrounds are the total 
SST changes during 1970-2006 for the nine regions. The 
green background behind each regional SST change indicates 
that the change is considered favorable for increased TC 
activity. 


2. Shear 

Equation 4 shows a surprisingly good least squares fit 
for shear as a function of time and location (as evidenced 
by R-squared value of 0.4862). 

Shear (m/s) = -16.04 + 0.0122year (4) 

where the standard errors for the intercept and the year 
coefficient are 9.02 and 0.0045, respectively. 
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This equation shows that shear, like SST, is dependent 
on year (our proxy for global warming) , and the 
contribution of year is about a 0.4 m/s increase over the 
past 36 years. The coefficient for our year term, +0.0122 
is an estimate of the relationship between year and shear, 
and has units of (m/s)/year. As year is our proxy for 
global warming, equation (4) indicates that there is a 
positive relationship between global warming and WNP shear. 
Note too, that again the coefficient scales the equation 
and keeps the units consistent, so the size of the 
coefficient is not indicative of its importance. 

Figure 14 shows that the shear is near zero or 
negative, particularly in the tropics. Moreover, the 
increase in shear is most readily seen in the lowest 
latitudes of our region of concern. Here "increasing 
shear" might be misleading because the shear increase means 
shear is tending to be less negative, but the magnitude of 
shear (which is what is important to TC genesis and 
sustainment), is actually decreasing. 

As we saw with SST, shear is changing with time, and 
its change is clearly regionally dependent. This supports 
our contention that analysis of the relationships between 
TC activity and the LSEFs needs to be done with data on the 
local level e.g., our weekly five by five blocks. 
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Figure 13. 
only as a 
increase in 


Least squares fit of shear data using year 
dependant variable. Note the very slight 
average shear with time. 
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Figure 14. Shear vs. year conditioned on latitude and 

longitude. The change in shear (1970-2006) is shown for 
each region with green background indicating the change is 
considered favorable for increased TC activity and red 
considered unfavorable. 


3. Vorticity 

Since absolute vorticity is in part due to Coriolis 
effects (a function of latitude), absolute vorticity lends 
itself nicely to using least squares. The equation for 
absolute vorticity is: 


Vorticity (s ^*10^) = -1.100 + 0.0004year + O.OOlVweek 

+ 0.2184maxLat (5) 

where maxLat is the upper latitude of the 5 degree by 5 

degree blocks (all data, in its geographic and weekly 
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a maximum 


blocks is identified with a year, a week, 
latitude and a maximum longitude). The standard errors for 
the intercept and each of the coefficients in order are: 
0.4232, 0.0002, 0.0002, and 0.0002. The value for adjusted 

is 0.90. Analysis of variance (ANOVA) indicates that 
latitude is indeed an important predictor for absolute 
vorticity (p-value less than 0.00001) . Year, in the 
presence of week and latitude according to ANOVA is also a 
significant predictor (though barely: p-value = 0.048) 

whose contribution over the past 35 years has grown a 
negligible 0.014 s~^. Figure 13 also demonstrates that 

absolute vorticity is predominantly a function of latitude 
only, and not of year, our measure of global warming. 
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Figure 15. Absolute vorticity vs. year, conditioned on 
latitude and longitude. Note the predictor variable that 
dominates the response is latitude. 


4. Mean Upward Vertical Velocity 

Using a least squares fit of vertical velocity vs. 

year we found that there is a negligible change in vertical 
velocity over time, our proxy for global warming (see 

Figure 16). 

Figure 17 also shows a lack of dependence on year in 

our plot conditioned on latitude and longitude. Note also, 
that vertical velocity changes in the transition from the 
Reanalysis I field (1970-1979) to Reanalysis II (1980- 

2006); this change is most obvious in the positive 
extremes. 
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Figure 16. Least squares fit of vertical velocity vs. 

year. 
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Figure 17. Box plot of vertical 

conditioned on latitude and longitude. 


velocity vs 


year 


5. Relative Humidity 

When examining relative humidity using a simple least 
squares fitted line, we found that there was a significant 
positive trend over time (Figure 18). However, examination 
of Figures 19 and 20 indicates that the upward trend is 
artificial and due to differences between the Reanalysis I 
and II data sets. 
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Figure 18. Least squares regression line based on both 

Reanalysis I and II data showing an increase in relative 
humidity over time. 

Relative humidity, using only data from the reanalysis 
II data set, was analyzed with the least squares fitted 
line shown in Figure 21. We see that when analyzing one 
consistent data set, relative humidity, instead of sizably 
increasing over time over the WNP, decreased very slightly 
(a total of about 0.37 percent (compared to about 70 
percent) over the past 25 years. Figure 22 indicates a 
seeming lack of dependence of relative humidity on latitude 
or longitude. 
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Figure 19. Box plot of relative humidity showing 
discontinuity at the 1979-1980 transition between the 
Reanalysis I and II data sets. 
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Figure 20. Box plots of relative humidity conditioned 
on latitude and longitude. The discontinuity between data 
sets at 1980 is again readily apparent. 
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Figure 21. Plot of least squares 

II relative humidity data only. 


fit using Reanalysis 
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Figure 22. Relative humidity box plots conditioned over 

latitude and longitude. 


In summary, we have found that of the LSEFs, SST and 
shear seem to be in part a function of year, our proxy of 
global warming, and that in general they have become more 
favorable for TC activity. Thus, SST and shear in the WNP 
may have been altered by global warming in a way that would 
favor increased TC activity. The other LSEFs do not appear 
to be functions of year, and thus: (a) their effects on TC 
activity have apparently not changed with time; and (b) 
they do not appear to have been altered by global warming. 
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B. LONG TERM TRENDS IN TROPICAL CYCLONE ACTIVITY 

1. Tropical Cyclone Frequency 

Figure 23 shows the least squares fitted line of TC 
frequency in the WNP during 1970 - 2005. Clearly there is 
an increasing trend with time, with the equation for this 
line as follows: 

TCCount = -273.82 + 0.1529*year (6) 

where year is considered significant as a predictor, having 
a p-value of 0.048. The standard errors for the intercept 
and the year term are 148.27 and 0.0746 respectively. 

Equation 6 indicates that the average number of TCs 
has increased since 1970 by over five per year, again over 
the entire WNP. for the regression is exceptionally low 
(0.11) . This low value is an indication that TC frequency 
is dependent on more than just year, but also other 
independent variables, such as the LSEFs. 
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Figure 23. Least squares linear regression of TC 

formation for 1970-2005 in the Western North Pacific. 


While there does appear to be some nonlinearity in the 
expected annual TC formations as a function of time, the 


variability 

of 

annual 

TC 

formations in 

Figure 

23 appears 

relatively 

constant 

over 

time. The 

normal 

quantile- 

quantile plot 

of the 

residuals of the 

regression fit in 


Figure 24 supports that the underlying regression 
assumptions are reasonable for this data fit. 
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Quantiles of Standard Normal 

Figure 24. Normal probability plot showing a normal 

distribution of the residuals. 

So far, we have seen that of the LSEFs, SST and shear 
seem to be in part a function of year, our proxy for global 
warming, and that for the most part they are tending in a 
direction which favors TC activity. We now also have 
Equation 6, which on a very broad analysis scale indicates 
that TC formations are also a positive function of our 
global warming proxy, year. 

2. ACE 

Eigure 25 demonstrates a clear upward trend in ACE 
over time in the Western North Pacific over the 1970 
through 2005 timeframe. Using linear regression we get the 
following equation: 
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ACE (kts^) = -52027255 + 27695*year (7) 

Equation 7 indicates that over the past 36 years, 
annual average ACE in the WNP has increased on average by 
almost one million knots^. Eigures 25 and 26 suggest that 
the linear regression modeling assumptions of 
homoschedasticity and normal errors are reasonable. 

As with the results for TC frequency, the value of R- 
squared is low (0.083), which is again due the fact that in 
this model the only independent variable in Equation (7) is 
year, and if we wanted to predict ACE with accuracy (rather 
than merely demonstrate the effect of our global warming 
proxy on ACE) we would need to include other independent 
variables such as the LSEEs. 

There is some slight evidence (p-value = 0.0883)) that 
the coefficient for year estimated in Equation (7) is not 
zero. Had the p-value been 0.05 or less, we would have 
claimed that ACE is a function of year, our proxy for 
global warming. Had the p-value been greater than 10% we 
would have stated that ACE does not appear to be a function 
of year. As it stands, we can only state that there may be 
a relationship, and try to establish what that relationship 
is when we examine ACE as a function of the LSEEs later in 
this study. 
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Annual Ace 



Figure 25. ACE (kt^) vs. year for the WNP, 

trend included. 


with 


linear 
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Quantiles of Standard Normal 

Figure 26. Normal probability plot of the ACE residuals 

showing the residuals are indeed normally distributed. 

The ACE results are intriguing. We expect that if TC 
formations increase, there might be a corresponding 
increase in ACE. However, Eigure 27, which is a plot of 
ACE per storm for each year, shows an increase in ACE over 
time even though this average should show a trend line that 
hovers around zero if the increase in ACE was due solely to 
an increased number of storms. This is consistent with 
other studies^? which suggest that the intensity of 
individual TCs has been increasing over time. 

We have shown that of the LSEEs, SST and shear seem to 
be in part a function of year, our proxy of global warming, 

P. J. Webster, et al.. Science, 309 , 1844 (2005). 
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and that for the most part they are tending in a direction 


which favors TC activity. We also have Equation 4, which 
on a broad scale indicates that TC formations are also a 
positive function of our global warming proxy, year, and 
Equation 7, which on a broad scale indicates that ACE is 
increasing. Eurthermore, it appears that ACE is increasing 
more than can be accounted for by the increase in the 
number of TCs. 

These results suggest (but do not prove) that: (1) 
global warming during 1970-2005 may have affected two of 
the LSEEs in the WNP, SST and shear, and (2) global warming 
may have contributed to an increase in TC frequency and ACE 
during 1970-2005. In addition, the changes in SST and 
shear are consistent with the TC frequency and ACE changes. 


indicating 

(but 

not 

proving) that global warming 

induced 

changes in 

SST 

and 

shear may 

have 

contributed 

to the 

increase in 

TC frequency and ACE 

during 

this period. 



53 



150000 - 


100000 - 
Ol 

o 

< 

> 

< 


50000 - 


0 -I— 
1960 


Figure 27. Average ACE per storm vs. year showing a 

small but pronounced upward trend. 

C. MODELLING OF TROPICAL CYCLONE ACTIVITY IN THE WESTERN 
NORTH PACIFIC USING LARGE SCALE ENVIRONMENTAL FACTORS 


1. Tropical Cyclone Frequency 

Using Equation 7, we have shown a relationship between 


year and 

the 

number 

of 

TCs 

formed. 

Based on 

using 

the 

total number 

of TCs 

per 

year 

in the 

WNP, 1970 

-2005, 

and 

ignoring 

the 

LSEEs. 

We 

have 

also shown that there 

is a 


relationship between global warming and two of the LSEEs. 
We now set about modeling TC formation where the 
independent variables are the LSEEs, the factors known to 
influence TC formations. 
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The variable used for TC formation is binary, that is 
either a TC formed in our weekly five by five degree block, 
or a TC did not form. Logistic regression was used as our 
modeling tool to estimate the probability of TC formation p 
as a function of the independent variables. We started 
simply, using SST as the only independent variable. 

One problem that immediately became apparent was that 
the model is significantly underdispersed. Underdispersion 
is indicated by the fact that the residual deviance should 
be approximately equal to its number of degrees of 
freedom^O; our results show a residual deviance of 10652.02 
on 144934 degrees of freedom. This underdispersion 
indicates that the variability of the response variable is 
much less than p(l-p), the variance of the Bernoulli 
distribution assumed in logistic regression. 

The cause of underdispersion might be explained 
because in logistic regression observations are assumed to 
be independent, whereas these may not always be. For 
example TC formation in a region means that some relatively 
slowly changing and wide spread influences on TC formation 
(such as SST) may favor the formation of other TCs. 
Conversely the formation and growth of a TC may alter 
conditions (such as shear or SST) so as to decrease the 
likelihood of TC formation in the immediate future and in 
the immediate proximity. Our conjecture is that this 
altering of conditions preventing other TCs from forming 
would act only through the LSEFs, and minimize this 
dependency problem. 


D. Collett, 1991, Modelling Binary Data. Chap. 5, Overdispersion. 
Chapman and Hall, 2-6 Boundary Row, London SEl 8HN, UK, 188. 
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with this underdispersion, caution must be taken when 
interpreting inference results based on the logistic 
regression fit to this data. However, the logistic 

regression fit to the data can still be used to estimate 
the probability of TC formation. Whether the model fits 
the data or not will be assessed by comparing how well 
estimated probabilities of TC formation compare to actual 
formations of an independent set of data. 


In 

a logistic 

regression model, the 

probability 

of 

success, 

P (in 

this 

case "success" is the 

formation of 

a 

TC) is 

linked 

to k 

independent variables 

Xi, X2, . . . Xk 

as 


follows: 

logit(;?) = p^ + + pj^x^ +... + p^x^ (8) 

where \ogii{p) , which is log (p/(1-p) ) , is the log-odds of the 
success probability^!, and the are coefficients of 

the linear combination of the k different independent 
variables. 

In its simplest version (the null model), the success 
probability is a function only of a constant term, . For 
progressively more complex models, the success probability 
is a function not only of the term, but of one or more 

independent variables and their corresponding coefficients. 

Using a combination of variable selection techniques 
and diagnostic plots gives the model fit summarized in 
Table 3. 


D. Collett, 1991, Modelling Binary Data. Chap. 5, Overdispersion. 
Chapman and Hall, 2-6 Boundary Row, London SEl 8HN, UK, 188. 
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Term 

Coefficient 

Standard Error 

Intercept 

-22.3677 

2.3205 

SSI 

0.0066 

0.0007 

Coriolis 

1.3588 

0.1954 

AbsVort 

0.0407 

0.0254 

VertVel 

0.1436 

0.0095 

RelHum 

-0.0276 

0.0097 

Shear 

-0.0264 

0.0060 

Shear^ 

-0.0031 

0.0004 


Table 3. Estimated coefficients and standard errors for 

model to predict the formation of a TC in a weekly 
five by five degree block. 

In the course of model development and variable 
selection we learned that our model was very sensitive to 
the changes between the two NCEP reanalysis data sets. As 
the majority of our data was from the better Reanalysis II 
data set, we used only the Reanalysis II data set (1980- 
2006) in the development of the model. At the same time, 
we reduced the data set, retaining only the weekly five by 
five degree blocks with TC formations or ACE, as well as a 
randomly selected twenty percent of the remaining TC 
activity free blocks. This was done to reduce computer 
runtime and memory issues, and it also reduced, but did not 
eliminate our underdispersion problem. 

Table 3 shows that we have included a Coriolis term in 

Coriolis effects are an implicit large scale 
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our model. 




factor that is necessary for TC formation^^^ with values 
near zero (e.g., within a few degrees of the equator) being 
unfavorable for TC activity. The Coriolis effect is a 
strong function of latitude and is unlikely to change 
significantly with global warming. However, when running 
the above model without a Coriolis term, the result was 
significant overprediction at and close to the equator. 
Thus, we included a Coriolis term in the model. 

Care must be taken when interpreting the magnitude of 
the estimated coefficients given in Table 3. These 
coefficients represent the partial effect of each predictor 
in the presence of the rest. This is more understandable 
when you examine Figure 28. There are noticeable, though 
not terribly strong relationships between mid-level 
relative humidity and in particular, SST and vertical 
velocity. These relationships obscure the contribution to 
TC formation by relative humidity, because, if you have 
information on SST for example, you also have some 
information on relative humidity. 


W. M. Gray, 1975: Tropical Cyclone Genesis in the Western North 
Pacific, Naval Air Systems Command, Washington DC, 20361. 
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Figure 28. Scatter matrix of time, main factors, CO 2 

concentration, and temperature anomaly. 


Table 3 also indicates that among the independent 
variables, a quadratic term of shear is also included, and 
is significant. Polynomial terms are necessary to model 
the log-odds as a non-linear function of shear. While shear 
ranges from approximately -40 m/s to 70 m/s, it is most 
conducive to TC activity at near zero values. A graph 
showing the contribution of the shear terms to the 
estimated log-odds is seen in Figure 29. The slight left 
skew of the shear contribution may be due to some physical 
phenomenon, but most likely is due to the fact that most 
TCs form in or near the tropics where the average shear is 
negative. It is no surprise that the model with the best 


59 









































fit demonstrates that the contribution of shear to the 
probability of TC formation is zero when shear is near 
zero, and becomes sharply negative as the magnitude of 
shear increases. 



Figure 29. Contribution to TC formation probability 
model by the main effect shear. 

The one obvious remaining issue is that the model of 
Table 3 indicates a negative relationship between relative 
humidity and TC formation, the opposite of what was 
expected. However, Figure 30 helps explains this issue. 
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Figure 30. Scatter matrix of relative humidity and 
vertical velocity. 

Figure 30 shows there is a strong positive 
relationship or multi-colinearity between relative humidity 
and vertical velocity. This multi-colinearity, combined 
with the knowledge that high relative humidities are 
required for TC formation and that vertical velocity has an 
upper bound on its contribution to TC formation, causes the 
coefficient for our relative humidity term to be weakly 
negative, even though the marginal contribution of relative 
humidity to TC formation is most definitely positive. The 
residual plots in the next few sections help shed light on 
this behavior. 
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The end result of Table 3 is that if we want to 


characterize the contribution of each of the LSEFs to TC 

formation, then an increase in SST, absolute vorticity, 

and/or vertical velocity - or a reduction in the 

magnitude of shear-will result in an increase in TC 

formation probability in the presence of the other factors. 

A more intuitive feel for the relationship between TC 
formation and the LSEFs can be obtained by examining the 
partial residual plots for each of the factors. Partial 

residual plots are used as (among other things) diagnostics 
to check the form of the independent variables. In this 
case each of the relationships between the response and an 
independent variable appears linear (except for having to 
account for crossing through zero with shear as explained 
earlier) , and so no further transformations of the data 
were required. 

a. The Effect of SST on TC Formation 

Figure 31 shows the partial residual plot for SST 
against TC formation. To aid the reader who is unfamiliar 
with partial residual plots, in Figure 31, the lower "line" 
of data represents data points for which a TC did not form 
for each weekly five by five degree block (at the specified 
SST (°C*100)) . The upper line is comprised of all the 
weekly blocks at the specified temperature where a TC did 
form. The vertical axis is a measure of the partial 
effects of SST on the log-odds that a TC will form. 

There are approximately 1000 TC formations 
represented in Figure 31. Notice that perhaps one TC out 
of 1000 formed in water for which the temperature was less 
than 26 °C. Clearly, as you examine the figure from left to 
right, examining the density of the data points, it is 
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obvious that the warmer the water is, the more likely it is 
that a TC will form. The residual plot of Figure 31 allows 


us to draw some important conclusions: 


1. The LSEF identified in Chapter I, sea surface 

temperatures (SSTs) exceeding 26 seems to 

be confirmed by our residual plot. 

2. Increasing SST increases the partial effect of 
SST on the log-odds that a TC will form. 

3. We have seen earlier that global warming appears 
to have a positive relationship with SST. Now 
we observe that increasing SSTs appear to have a 
positive relationship with TC formation. 



SST Deg C*100 


Figure 31. Partial residual plot for SST as used in our 

logistic regression model. 


W. M. Gray, 1975: Tropical Cyclone Genesis in the Western North 
Pacific, Naval Air Systems Command, Washington DC, 20361. 
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b. The Effect of Shear on TC Formation 

Figure 32 is the partial residual plot for our 
shear term. Again, a clear pattern can be observed. Out 
of over one thousand tropical cyclones examined, fewer than 
ten formed in conditions where the magnitude of the shear 
exceeded 20 m/s. The vast majority of TCs which form do so 
where the magnitude of the shear is less than 20 m/s, with 
the greatest concentration occurring exactly where one 
would expect, where shear is near zero. 

The residual plot of Figure 32 allows us to draw some 
important conclusions: 

1. The LSEF identified in Chapter I, weak vertical 
shear of the horizontal winds^^, seems to be 
confirmed by our residual plot. 

2. The occurrence of vertical shear with a magnitude 
greater than 20 m/s means it is very unlikely a TC 
will form. 

3. We have shown in prior sections that for the WNP, 

and based on 1970-2005 data: (a) global warming 
appears to have a positive relationship with TC 
frequency; (b) global warming appears to have a 
positive relationship with shear; (c) global warming 
appears to have a positive relationship with shear 
magnitude. These results, combined with the 

positive relationship between low shear magnitude 
and TC frequency (Figure 32) indicates that the 
apparent impacts of global warming on TC frequency 
in the WNP have occurred via the effects of global 
warming on shear. 


W. M. Gray, 1975: Tropical Cyclone Genesis in the Western North 
Pacific, Naval Air Systems Command, Washington DC, 20361. 
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Figure 32. Partial residual plot for shear as used in our 
logistic regression model. 


c. The Effect of Relative Humidity on TC 
Formation 

We have hypothesized that relative humidity plays 
a role in TC formation. The partial residual plot of 
Figure 33 bears this out, as almost all TCs that have 
formed in the Western North Pacific since 1970 did so when 
the relative humidity was 55% or greater. Certainly Figure 
33 demonstrates that high relative humidity is positively 
associated with TC formation. The greatest concentration of 
TC formations occurs at what appears to also be the 
greatest concentration of relative humidities 
(approximately 75 to 80%). This suggests that relative 


65 





humidity plays a threshold role; that is, once relative 
humidity exceeds a certain limit, perhaps 55%, the relative 
humidity is favorable to TC formation, but does not further 
increase the formation probability as humidity continues to 
climb. Clearly Figure 33 confirms the LSEF: high mid-level 
humidity^^. 



Figure 33. Partial residual plot for relative humidity 

as used in our logistic regression model. 


W. M. Gray, 1975: Tropical Cyclone Genesis in the Western North 
Pacific, Naval Air Systems Command, Washington DC, 20361. 
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d. The Effect of Vertical Velocity on TC 
Formation 

For this discussion of vertical velocities, 
recall from Chapter II that we multiplied the NCEP 
reanalysis vertical velocities (in (Pa/s) by -1 prior to 
analyzing them. Thus, in our results, positive vertical 
velocities correspond to upward motion of the atmosphere. 

Figure 34 shows a clear relationship between 
vertical velocity and TC formation. Not one TC formed in 
conditions where the vertical velocity was below -5 Pa/s 
and few formed above approximately 25 Pa/s. Thus, the LSEF 
mean upward velocity^® is borne out by Eigure 34. 


W. M. Gray, 1975: Tropical Cyclone Genesis in the Western North 
Pacific, Naval Air Systems Command, Washington DC, 20361. 
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Figure 34. Partial residual plot for vertical velocity 

as used in our logistic regression model. 

e. The Effect of Absolute Vorticity on TC 
Formation 

We hypothesized that absolute vorticity is 
significant as a factor which influences TC activity. 
Figure 35 helps confirm the LSEF, large positive absolute 
vorticity at low levels'^'^, as it shows that virtually all 
TCs have formed when absolute vorticity was greater than 
zero. Of note, it appears that the concentration of TC 
formations decreases sharply at approximately 9 s~^10^. 

Absolute vorticity has a very strong relationship to 
latitude, and most likely the decline in the concentration 

W. M. Gray, 1975: Tropical Cyclone Genesis in the Western North 
Pacific, Naval Air Systems Command, Washington DC, 20361. 
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of TC formations with higher absolute vorticity is actually 
due to the fact that higher vorticity is associated with 
higher latitude and thus, lower SST and higher shear. 
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It is highly unusual for a TC to form when 
absolute vorticity (in a weekly 5 by 5 
block) is negative or greater than 9 s-’1CP 
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higher latitudes 
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Figure 35. Partial residual plot for absolute vorticity 

as used in our logistic regression model. 


f. Model Validation 

Thus, the model of Table 3 leads to the equation 
below which is used to estimate the probability of TC 
formation: 

logit(p)= -22.367 + 0.0066SST + 1.359Coriolis 
+ 0.0407AbsVort + 0.1436VertVel 
-0.027GRelHum - 0.0264Shear 
- 0.0031Shear2 (9) 
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All the terms in Equation 9 are physically 
reasonable, given the explanation of the negative sign for 
the relative humidity term presented previously. However, 
as stated earlier when we first discussed the 
underdispersion of our model, we must demonstrate model 
validity by how well it estimates the probability of TC 
formation. We constructed contour plots of estimated 
probability for each of our weeks included in the data base 
for the years 1980 through 2005. If the model is valid, it 
ought to be able to accurately map out regions of the WNP 
with relatively high TC formation probabilities. By 
plotting actual TC formation locations for each week on the 
probability contour plot for that week, we can assess the 
accuracy of the model probability maps. Figures 36-38 show 
examples of model weekly TC formation probability maps 
along with actual TC formation locations. Note the 
resemblance of the higher probability regions to regions 
known to be generally favorable for TC formation (e.g., the 
monsoon trough at about 5-15N). 
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Figure 36. Model weekly probability contours for TC 

formation for week 43 of 1997. Red dots mark formation 
locations for actual TCs. 
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Figure 38. Model weekly probability contours for TC 
formation for week 50 of 1997. Red dots mark formation 
locations for actual TCs. Note the low probabilities 
compared to those in Figures 36-37. Note also that no TCs 
formed during this week. 

While the Figures 36-38 indicate that the model 
is valid, the results do not constitute a true check of 
model validity: the model was built on data from 1980 
through 2005, and thus should have good performance over 
those years. A more definitive check uses the estimated 
model to predict TC formation for years not included in the 
data set (similar to hindcasting in atmospheric 
forecasting). Even though we had excluded the 1970 through 
1979 data from our model (because of differences between 
the data sets), it stood to reason that our model should be 
able to at least roughly estimate TC formation 
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probabilities during those years. Figures 39-42 show the 
performance of our model in a sampling of weeks in 1974. 
Clearly our model works on data outside of the data set 
used for model development, thus substantiating the 
validity of the model. 




Figure 39. Model weekly probability contours for TC 

formation for week 43 of 1974. Red dots mark formation 
locations for actual TCs. TC data for 1970-1979 was not 

used in the model development. Thus, the results in this 

figure and in figures 39-42 indicate that there is a good 
fit between our model and actual TC formation in the 

Western North Pacific. 
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Figure 40. Model weekly probability contours for TC 

formation for week 42 of 1974. Red dots mark formation 
locations for actual TCs. TC data for 1970-1979 was not 

used in the model development. Thus, the results in this 

figure and in figures 39-42 indicate that there is a good 
fit between our model and actual TC formation in the 

Western North Pacific. 
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Figure 41. Model weekly probability contours for TC 

formation for week 39 of 1974. Red dots mark formation 
locations for actual TCs. TC data for 1970-1979 was not 

used in the model development. Thus, the results in this 

figure and in figures 39-42 indicate that there is a good 
fit between our model and actual TC formation in the 

Western North Pacific. 
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Figure 42. Model weekly probability contours for TC 
formation for week 20 of 1974. Red dots mark formation 
locations for actual TCs. TC data for 1970-1979 was not 
used in the model development. Thus, the results in this 
figure and in figures 39-42 indicate that there is a good 
fit between our model and actual TC formation in the 
Western North Pacific. Note the low probabilities compared 
to Figures 39-41 and the absence of any TC formations. 


The probability maps in Figures 36-42 constitute 

some of the better matches between model probabilities and 

actual TC formations that we obtained. However, the 

overall results do demonstrate sufficient goodness of fit 

that we can claim that the relationships indicated by our 

model are valid. Appendix A and Appendix B have all the 

weekly probability contour plots from 2004 and 1979 

respectively so that the reader can estimate for him or 

herself how useful a predictive tool the above model is. 

The model appears to have potential as a tool for 
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predicting the number and location of TC formations (note 
though that the development of such a tool was not a goal 
of our study) . However, for use as a predictive tool, the 
model needs further development, and more quantitative and 
formal assessments of its predictive skill, as discussed in 
Chapter IV. 

g. Sensitivity of Model Probabilities to Each 
LSEF 

Once we had a model that showed the relationship 
between TC formation probability and the LSEFs, we decided 
to investigate the sensitivity of the model probability to 
each LSEF, as a way of estimating the relative contribution 
of each LSEF to the probability of TC formation. To 
demonstrate this relative contribution, we selected what 
appear to be representative values for each of the LSEFs at 
the times and locations at which the model indicated high 
probabilities of TC formation. We then altered the LSEFs 
by moderate amounts (ten percent of the typical range of 
the LSEF over which TC formation was likely to occur) while 
holding the other factors constant, to see the effect on 
formation probability. The results are (shown in Table 4) 
indicate that the model's TC formation probabilities are 
most sensitive to changes in vertical velocity and SST, and 
least sensitive to changes in absolute vorticity and 
relative humidity. 
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Factor 

Normal 

Base 

10% 

Base formation 

Probability 

Percentage 


range 

value 

Change 

probability 

change 

change 

SST 

2600- 

2800 

+ 40 

0.07614 

0.0208 

+27.28 

(deg C*100) 

3000 






Abs Vort 

2-10 

6 

+ 0.8 

0.07614 

0.0023 

+ 3.05 

(s^^*ioh 







Vert Vel 

0-20 

5 

+ 2 

0.07614 

0.0228 

+29.97 

(Pa/s) 







RelHum (%) 

60-90 

70 

+ 3 

0.07614 

-0.0056 

-7.38 

Shear (kts) 

-20-20 

0 

+ 4 

0.07614 

-0.0102 

-13.43 


Table 4. Relative contribution to formation probability 

given a 10% change in each of the LSEFs while holding 
the other LSEFs constant. 


h. 

Relationship 

Formation 

Between Global 

Warming and 

TC 

As 

a reminder to 

the reader, we 

have shown 

the 


following: 

1. Global warming appears to have a positive 
relationship with two of the LSEFs, SST and shear. 

2. The changes in SST and shear have been in a 
direction considered favorable for TC formation and 
intensification. 

3. A model for TC formation probability confirms that 
SST and shear are necessary independent variables 
for modeling TC formation probability. 

4. The model indicates as SST increases or as shear 

tends toward zero, TC formation probability 

increases. 

5. These results indicate that global warming has a 
positive relationship with TC formation probability. 

6. Thus it appears that global warming induced changes 
in WNP SST and shear during 1970-2005 have 
contributed to an increase in TC formations during 
that period 

2. ACE 

To model the effect of the LSEFs on ACE, we were 
forced to take a slightly different approach. Examination 
of residual plots of linear regression fits and partial 
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residual plots such as Figure 43 for SST, reveal that the 
variability of ACE increases with its expected value. Thus 
the modeling assumption of equal variance in multiple 
linear regression models is violated. Note that the ACE we 
analyzed occurred in a weekly five by five degree block, 
and it is not the total ACE for a year or the ACE for a 
single TC. However, just as partial residual plots for TC 
formation were informative, they are informative for ACE as 
well, as shown in the following section. 

a. Developing a Predictive Model for ACE 
Although our data does not meet the assumptions 
necessary to use analysis of variance, we chose to fit with 
model development and then check the resulting fit against 
an independent set of data to demonstrate the validity of 
the model, as we did with the model for TC formation. The 
data was all from weekly five by five degree blocks which 
had either positive values for ACE, or were zero because a 
TC had formed but not yet gathered sufficient strength 
(maximum sustained winds of 18 m/s, or 35 kts) . This data 
set was composed of over 4500 data points. The model based 
on this data is shown in Table 5. 
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Term 

Coefficient 

Standard Error 

Intercept 

-92365.59 

9113.58 

SST 

26.52 

2.62 

AbsVort 

3005.25 

198.86 

VertVel 

58.83 

70.81 

RelHum 

211.99 

59.30 

Shear 

55.32 

42.36 

Shear^ 

-5.09 

1.49 


Table 5. Model coefficients and standard errors for ACE. 

Table 5 yields the following predictive equation 

for ACE: 

ACE (kt^) = -92365 + 26.52SST + 3005.25AbsVort 
+ 58.83VertVel + 211.99RelHum 
+ 55.32Shear -5.09Shear^ (10) 

The ACE that is predicted by the coefficients and 
independent variables shown in Table 5 and Equation 10, is 
the ACE that is predicted to be generated in a weekly five 
by five degree block, assuming a TC has formed in or is 
present in that block. Again the variability of ACE 
increases with its expected value as shown in the partial 
residual plots. This forces us again to check the model 
fit on an independent set of data to make sure the model is 
valid. Thus, as we did with the TC formation data, we have 
excluded the data prior to 1980. The ten years of data not 
used in the modeling proved to be very useful for 
determining goodness of fit. Notice also in Table 5 that 
the Coriolis term is no longer part of the model as it was 


81 




found to not be significant. This is not surprising, since 
this model describes ACE given that a TC (which is 
dependant on Coriolis acceleration) has already formed. 

We expected that the same independent variables 
that would be important for TC formation would be important 
for the generation of ACE. Certainly Table 5 bears this 
out, since the independent variables that are important for 
formation are also important for ACE (except for the 
Coriolis term). The major difference is that the strong, 
positive relationship between ACE and relative humidity is 
much more straightforward than the relationship we observed 
for TC formation. Also, exactly as expected, and consistent 
with our TC formation model, increasing SST, humidity, 
vorticity, and vertical velocity (up to our observed 
threshold limits) all tend to increase ACE, while any 
increase in the magnitude of shear tends to reduce ACE. 

We then proceeded to check our model that 
predicts ACE for a block given in which a TC had formed and 
was present, using our LSEE data as input for the model. 
We compared the annual actual ACE for each year to the 
annual predicted ACE for each year. The results are shown 
in Eigure 43. 
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Figure 43. Annual ACE values for 1970-2005. Actual (in 

blue) and modeled (in red) . 


Note that included in Figure 43 are actual (in 
blue) and predicted (in red) values for ACE for the years 
1970 through 1979. The data from these years was not 
included in the model fit. Moreover, as noted earlier in 
this chapter, there are some rather significant differences 
between the Reanalysis I and Reanalysis II data sets. 
Thus, it is encouraging that the predicted values of ACE 
match actual ACE rather closely during 1970-1979. 

Finally, Figure 43 seems to validate our model 
for ACE. There are clearly some differences between the 
model and actual values of ACE, but for the most part, both 
the magnitude and the trend (increasing or decreasing from 
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year to year) of the model seems to be accurate enough to 
state with confidence that the model works and that the 
relationships between the LSEFs and ACE hold true. Of 
course, additional quantitative assessment of the model 
skill would be useful. 

b. ACE Partial Residual Plots - SST 

Figure 44 shows the partial residual plot for ACE 
against SST. The straight horizontal line centered on zero 
represents all the weekly five by five degree blocks for 
which there was no TC activity (on the order of 140, 000 
data points) . All data above the horizontal line 

represents the value of ACE in all the weekly five by five 
degree blocks for which there was TC activity at its 
corresponding SST. ACE will be zero if a TC did not form 
(in some cases ACE is zero the week a TC formed if it did 
not gather enough strength to reach 35 kts prior to week's 
end) . 

Note that the ACE residual plot for SST has a 
different pattern than the corresponding figure for TC 
formation (compare Figures 44 and 31) . In particular, the 
non-zero ACE cases occur in SST ranges that are different 
from, but similar to, those for TC formations. This is 
because there is not a one for one correspondence between 
the weekly five by five blocks for TC formation and for 
ACE, since formations occur when the TC is below the 
intensity threshold for the calculation of ACE (maximum 
sustained winds of 18 m/s or 35 kts; see Chapter II) . 
Also, the non-zero ACE cases are spread over a wide range 
of the vertical axis. That is, the magnitude of ACE is a 
function of SST (as shown by Figure 44) . So not only do 
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higher SSTs make it more likely there will be a TC, they 
make it more likely that the strength of the TC will be 
greater too. 
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Figure 44. ACE partial residual plot for SST. 

c. Shear Partial Residual Plot for ACE 

For shear and the remainder of the LSEFs the 
magnitude of ACE can be quite large over a relatively broad 
range of LSEE values, so long as these values are at least 
near the thresholds identified in the TC formation residual 
analyses (see figures 45 through 48) . For example, as 
shown in Figure 45, as long as the value for shear is 
between -20 and 20, the value of ACE can be quite large or 
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it can be quite small (in fact it is skewed toward small) . 
But large values of ACE can and do occur anywhere within 


that range of shear. As we saw with TC formation, ACE 
peaks near a shear magnitude of zero. 



-^10 -20 0 20 ^10 60 8 

Shear 


Eigure 45. ACE partial residual plot for shear. 

d. Relative Humidity Partial Residual Plot for 
ACE 

Similar to the results for TC formation, a 
relative humidity of about 55-90% appears to be favorable 
for the generation of ACE (see Eigure 46). 
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Figure 46. ACE partial residual plot for relative 

humidity. 

e. Vertical Velocity Partial Residual Plot for 
ACE 

Figure 47 shows that a zero or positive value of 
vertical velocity is normally necessary to generate ACE. 
It is intriguing again that, as with TC formation, it seems 
that there can be too much vertical velocity, as the ACE 
residuals taper off sharply for vertical velocities greater 
than 22 Pa/s. 
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Figure 47. ACE partial residual plot for vertical 

velocity. 

f. Absolute Vorticity Partial Residual Plot for 
ACE 

As we saw with TC formation, the peak value for 
absolute vorticity for which most ACE occurs at is about 5- 
6 s~^*10^ (see Eigure 48) . ACE probably decreases beyond 

this peak as explained earlier: a high absolute vorticity 
most likely means the data is from a higher latitude, where 
there are lower SSTs and higher shear magnitudes, both 
unfavorable to TC formation and intensification. 
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Figure 48. ACE partial residual plot for absolute 

vorticity. 

g. Sensitivity of Model ACE to the LSEFs 

We investigated the sensitivity of the model 
predicted ACE to the LSEEs, as a way of estimating the 
relative contribution of each LSEE to the generation of ACE 
in a weekly five by five degree block. We selected the 
same representative values for each of the LSEEs as we did 
for TC formation (the range and peaks for the LSEEs are 
consistent for both TC formation and ACE) , and altered them 
by reasonable amounts (again, 10%) to see the effect on ACE 
generation. The results are shown in Table 6, with the two 
largest contributors being (in order) absolute vorticity 
and SST. Note that shear is a very slight positive 
contributor to ACE if the magnitude of shear is near zero. 
If the magnitude of shear grows significantly from zero, 

shear becomes an increasingly negative contributor to ACE. 
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Table 6. Relative contribution to ACE given a 10% change 

in each of the LSEFs while holding the other LSEEs 
constant. 

h. Relationship Between Global Warming and ACE 

In summary, we have shown the following: 

1. Global warming appears to have a positive 
relationship with two of the LSEEs, SST and 
shear. 

2. The changes in SST and shear are in a direction 
favorable for the generation of ACE. 

3. A model of ACE based on the LSEEs confirms that 
SST and shear are necessary independent variables 
for modeling ACE. 

4. The model indicates as SST increases or as shear 
tends toward zero, ACE tends to increase. 

5. Therefore, it appears that global warming has a 
positive relationship with ACE. 

6. Thus, it appears that global warming induced 
changes in WNP SST and shear during 1970-2005 
have contributed to an increase in ACE during 
that period. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


A. TRENDS IN THE LARGE SCALE ENVIRONMENTAL FACTORS 

We showed in Chapter III that WNP SSTs have increased 
during 1970-2006 (Equation 3) . We also saw that WNP shear 
has increased, but that the absolute value of shear has 
actually decreased, particularly in the tropics (Equation 
4). These assessments are on a large scale (approximately) 
ten by forty degrees, and on a seasonal timeframe. But if 
the LSEEs really do influence TC activity, then two of the 
LSEEs have changed in directions that would tend to 
increase TC activity, while the others two have been 
essentially unchanged. Therefore, an increase in TC 
activity would be expected in the data set. In addition, 
one might expect as greenhouse gasses continue to 
accumulate in the atmosphere that the observed trends in 
the LSEEs would continue as well. 

Eor SST and shear, we observed that the changes 
occurring over 1970-2006 were consistent with, and 
apparently attributable to, global warming. However, these 
changes did not occur uniformly across the entire WNP 
(Eigures 12 and 14) . Instead, they changed by different 
magnitudes and in different directions depending on the 
region of the WNP. Thus, analyzing the LSEEs by averaging 
the LSEEs over large geographic regions or periods of time, 
as done in prior studies of the impacts of global warming 
on TC activity, smoothes out these sub-basin scale 
variations and potentially eliminates the ability to detect 
important effects of the LSEEs on TC activity. 


91 



B. TRENDS IN TROPICAL CYCLONE ACTIVITY 

The observed long term changes in the LSEFs suggests 
that TC formations and ACE should have also increased 
during 1970-2005. Indeed, that is the case as we see an 
increase in actual TC formations, and an increase in ACE as 
well (Figures 23 and 25) . It is important to realize that 
there is an increase in ACE that is independent of that due 
merely to the fact that the number of TCs is increasing. A 
key point is that the average number of TCs per year has 
increased during 1970-2005, and that ACE per TC on average 
has increased as well (Figure 27). 

We have also identified a series of relationships 
between global warming measures, the LSEFs, and TC activity 
measures. We have seen a strong linear relationship 
between atmospheric greenhouse gasses and GASTA, one of our 
measures of global warming (Figure 3). We have also 
determined that there is a relationship between global 
warming and the LSEFs (Figures 10, 13, 15, 17 and 22) : SST 
is increasing, and shear is increasing (corresponding to a 
decrease in the magnitude of shear in the tropics) . The 
remainder of the LSEFs have remained essentially unchanged. 
We have confirmed that the LSEFs (SST, shear, vertical 
velocity, relative humidity, and absolute vorticity) do 
influence TC activity, both formation and ACE (see Chapter 
ITT sections C.l and C.2) . This suggests that if Earth 
continues to warm, tropical cyclones in the WNP will become 
on average more frequent, and on average more powerful. 
However, the patterns of SST and shear changes in the WNP 
indicate that future increases in shear may eventually 
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cause tropical regions of weak negative shear to become 
regions of positive shear, thus potentially offsetting the 
effects of increasing SST. 

It has been shown that TC activity is indeed a 
function of all the LSEFs. To attempt to characterize TC 
activity with only one or two of the factors is like 
assessing the value of a house based only on its number of 
bathrooms. Likewise, assessing the impacts of global 
warming on TC activity using only one of the LSEFs (e.g., 
SST, as in almost all prior studies of the impacts of 
global warming) is likely to produce misleading results. 

It has been shown that as a factor, SST is unique in 
that it is the only factor for which an increase raises 
formation probability and ACE values without apparent bound 
(Figures 31 and 44). All other LSEFs appear to have limits 
beyond which TC activity begins to decrease. For example, 
there are values of absolute vorticity and vertical 
velocity beyond which TC activity tends to decrease 
(Figures 34, 35, 47, 48). Non-zero shear tends to lower TC 

formation probability and ACE (Figures 32 and 45). 

We determined that for a same percentage change in any 
of the LSEFs, the factor with the greatest effect on TC 
formations is vertical velocity, followed by SST (Table 4). 
Absolute vorticity comes in third though at only about 10% 
of the contribution that SST and vertical velocity make. 
Relative humidity has a slight negative contribution to 
formation, and as expected, any shear whatsoever reduces 
formation probability. 

We also determined that for a same percentage change 
in any of the LSEFs, the factor with the greatest effect on 
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ACE is absolute vorticity, followed by SST (Table 6). 
Relative humidity, a positive contributor to ACE comes in 
third. In addition, ACE is not nearly as sensitive to non¬ 
zero shear as TC formation is. 

The common belief is that the LSEEs are necessary but 
not sufficient to cause or sustain TC activity. However, 
the probability contour plots matched with actual TC 
formations, hints that the LSEEs, in particular, high SSTs, 
may be both necessary and sufficient for TC formation, at 
least in the WNP. This idea of course conflicts with the 
well established concept that TC formation requires some 
combination of both the necessary LSEEs, and a TC 
triggering process (e.g., an easterly wave, Elsberry 1995). 
To test the idea that the LSEEs may be both necessary and 
sufficient, more rigorous quantitative assessments of the 
TC formation model's skill need to be conducted. 

The results of Chapter III, sections C.l and C.2, show 
that the LSEE conditions favorable for TC formation and ACE 
in the WNP should be slightly altered for some of the 
LSEEs, leading to the following: 

• SST: 26 °C or greater (unchanged). 

• Vertical wind shear of the horizontal winds: 

weak, where "weak" is defined as a magnitude of 
20 m/s or less. 

• Absolute vorticity: zero or positive but less 

than 13 s~^*10^ (in the northern hemisphere). 

• Vertical velocity: zero or upward but not 

exceeding -30 Pa/s. 

• Relative humidity at mid levels: in excess of 

55%. 
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Note that none of these changes are radical departures 
from the LSEFs defined early in this study. However, our 
results allow us to quantify the LSEFs more carefully based 
on more extensive data sets and analysis methods than in 
prior studies. For example, our results allow such 
descriptions of TC activity favorable conditions as "large 
levels of absolute vorticity" and "weak vertical shear" to 
be more carefully quantified and to be bounded by 
observationally based thresholds. 

We found that predicting formation probabilities and 
ACE values was very sensitive to changes in the generation 
of the reanalysis fields, and as such. Reanalysis I data 
needed to be excluded from the model formulation, though it 
was of sufficient quality to be used in quality of fit 
tests following model development. Extending Reanalysis II 
data back in time to at least 1970 would be very helpful 
for future research. 

The TC formation probability contour plots show that 
TC formation is a function of "right here, right now" LSEF 
data. That is, to attempt to characterize TC formation 
with seasonal and/or basin wide averaged values 
significantly dilutes the contribution each LSEF makes 
towards TC formation probability. And, since two of the 
LSEFs appear to be affected by global warming, it is 
important to account for at least weekly and five by five 
degree variations in the LSEFs when trying to assess the 
impacts of global warming on TC activity in the WNP. 

It was never our intention to develop a model so 
accurate that it could be used for forecasting the number 
and locations of TC formations. However, the probability 
contour plots (Figures 36 through 42) raise the possibility 
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that statistically based forecasting of TC formation of the 
sort used in this study may indeed be possible, and with 
reasonable accuracy. 

C. LESSONS LEARNED 

Assumptions were made regarding what the measures of 
TC activity should be, and as many authors in the 
literature search used ACE, we used ACE as well for the 
sake of consistency. Using wind speed squared when wind 
speed is 35 kts or greater makes ACE a rather nonlinear 
measure of the strength of a TC. In addition to the 
previously noted differences between how wind speed is 
measured by organization, calculating ACE without the storm 
radius (a measure normally missing from the best track data 
set) makes ACE a less than ideal measure of TC intensity. 

In hindsight, both ACE and PDI^® (power dissipation 
index - a measure similar to ACE that uses the cube of wind 
velocity instead of the square) may not be the best of 

measures for what humanity cares about. ACE and PDI may 

quantify the intensity of a TC over its life, but what 
humans generally care more about is the power of a storm at 
landfall and for the next 12 to 72 hours. If a strong 
storm churns away at sea for weeks at a time but never 
makes landfall, ACE for the season will rise dramatically 
but with the exception of the hazards to aviators and 

mariners at sea, the impact of that TC on humanity will be 
negligible. 

Keep in mind also that TC damage on landfall 

(monetarily) is determined largely by population density, 
wealth, and planning^^, and not by the size or intensity of 

K. Emmanuel, Nature, 436 , 686-688 (2005). 

R. A. Pielke Jr, C. Landsea, M. Mayfield, J. Laver, 

Bull.Amer.Meteor.Soc 86 , 1571 (2005). 
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the TC. One possible improvement on ACE and PDI is 
integrated kinetic energy (IKE) developed by Powell and 
Reinhold^^. 

We think that the models developed in our study 
provided more accurate assessments of the impacts of global 
warming on TC activity than those of prior studies. This 
is in large part because we did not use basin wide and 
seasonal average values for the LSEEs and TC activity, but 
instead used weekly five by five blocks of data. However, 
we used monthly mean SSTs for our raw SST data (for the 
reasons given in Chapter II). This monthly mean constraint 
may have restricted what our models could do, especially 
since SST appears to be a very important LSEE. In future 
studies, it would be worth the effort to interpolate the 
monthly data into weekly averages. Another option 
available is the NCEP CISST data set (six hourly and 1 
degree resolution) . We had rejected this data set as an 
option because it only extended back to 1981 (Table 1) . 
But since our models ended up being based on Reanalysis II 
data back to 1980, this limit at 1981 now does not seem as 
severe as when we began exploring data sets. 

This research would have been far easier had the lead 
author taken courses in data mining and non-parametric 
statistics. Should others pursue a follow on topic from 
this one, these courses would be a worthwhile investment of 
their time. 

D. RECOMMENDATIONS FOR FUTURE STUDY 

This work just barely begins to scratch the surface of 
the wealth of information that can be gleaned from our data 


M. D. Powell, T. A. Reinhold, Bull .Amer.Meteor. Soc 88 , 513 (2007) . 
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set. Subjects for future study should include, but not be 
limited to, the following questions and topics: 

1. Can similar modeling be done in other regions 


such as the 

Atlantic? 

Would 

the same 

results be 

generated? 






2. How similar are 

the 

TC 

formation 

and ACE 

generation 

processes, 

as 

revealed 

by the 

regression 

models, 

for 

the 

different 

tropical 


basins ? 

3. Can similar modeling be done using the JMA and 
other TC data sets? Would the same results be 
generated? 

4. Use regression to show a relationship between 
JTWC data and JMA data (as well as other TC data 
bases). The differences between existing TC data 
sets are adding a lot of confusion as to what the 
long term patterns are in TC activity. 

5. We have seen that temporal changes in the LSEFs 

are region dependant. Can we quantify those 

changes? Can we predict those changes for other 
parts of the planet? 

6. Refine the models developed in this work. If at 
all possible, use a quality SST data set at a 
higher temporal resolution. 

7. Is it possible that the TC formation model or the 
ACE model could be used as an aid for medium and 
long range forecasting of TC formations and 
intensity? 
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APPENDIX A 


TC FORMATION PROBABILITY CONTOUR PLOTS 
FOR WEEKS 20-52 OF 2004 
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APPENDIX B. TC FORMATION PROBABILITY CONTOUR PLOTS 

FOR WEEKS 20-52 OF 1979 
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